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ABSTRACT 

The  present  work  is  concerned  with  the  nature  of  the  interior  solution 
and  the  influence  coefficients  of  shallow,  spherical,  thin,  elastic 
shells  (or  equivalently  a  shallow,  thin,  elastic,  paraboloidal  shell  of 
revolution)  which  are  homogeneous,  isotropic,  closed  at  the  apex, 
and  of  uniform  thickness.  The  investigation  is  carried  out  within  the 
framework  of  the  usual  shallow  shell  theory  for  small  displacements 
and  negligible  transverse -shear  deformations.  Exact  interior  solu¬ 
tions  are  obtained  for  shells  acted  upon  by  edge  loads  and  edge  mo¬ 
ments.  The  constants  of  integration  associated  with  these  interior 
solutions  are  expanded  asymptotically  in  inverse  powers  of  a  large 
parameter.  Retaining  only  the  leading  term  of  these  expansions  leads 
(in  most  cases)  to  known  approximate  results.  Explicit  expressions 
for  the  second-order  terms  are  obtained.  It  is  shown  that  these 
second-order  terms  play  a  significant  role  in  a  certain  class  of  prob¬ 
lems.  The  relative  importance  of  the  membrane  and  inextensional 
bending  stresses  in  the  interior  of  the  shell  is  discussed.  The  exact 
and  asymptotic  influence  coefficients  are  obtained.  The  interior 
stress  state  of  shells  subjected  to  polar  harmonic  axial  surface  loads 
is  also  investigated  by  the  same  procedure. 
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NOMENCLATURE 


Cartesian  coordinates  in  base  plane  (plane  tangent  to  apex 
of  undeformed  middle  surface  of  shell)  with  origin  at  apex 

z 

Distance  between  base  plane  and  a  point  on  undeformed 
middle  surface 

r 

Distance  between  a  point  on  undeformed  middle  surface  to 
axis  of  revolution 

e 

Polar  angle  in  base  plane  measured  from  positive  x-axis 

Circumferential  base  vector  of  undeformed  middle  surface 

Meridional  base  vector  of  undeformed  middle  surface 

n 

Unit  normal  of  undeformed  middle  surface,  positive  inward 

F 

Stress  function  from  which  in-plane  stress  resultants  are  derivable 

w 

Transverse  middle  surface  displacement 

u 

Meridional  middle  surface  displacement 

V 

Circumferential  middle  surface  displacement 

N  .  N  N 

iNr, 

In-plane  stress  resultants 

Mf’  ^0’  ’^re 

Moment  resultants 

Qr’  % 

Transverse-shear  resultants 

^r’ 

Rotations  of  middle  surface 

p 

n 

Normal  component  of  surface  load  intensity 

p 

r 

Meridional  component  of  surface  load  intensity 

Circumferential  component  of  surface  load  intensity 

P 

z 

Component  of  surface  load  intensity  in  direction  of  axis  of 
revolution,  i.e.,  the  axial  component 

D 

Bending  stiffness  factor 

1/A 

Stretching  stiffness  factor 

Poisson’s  ratio 

h 

Wall  thickness  of  shell 

E 

Young’s  modulus 

a 

Value  of  r  at  edge  of  shell 

(), 

_9() 

9r 

^\e 

_  9() 
ae 

v^O 

=  ()rr  +  l/r()r  +  l/r^),ee 

Vll 


Bending  stress 


Direct  or  membrane  stress 

Particular  solution  for  quantity  in  bracket 

Radius  of  spherical  middle  surface  or  twice  focal  length 
of  paraboloidal  surface  of  revolution 


MEMBRANE  AND  BENDING  STRESSES 
IN  SHALLOW,  SPHERICAL  SHELLS 


I.  INTRODUCTION 

One  of  the  characteristic  behaviors  of  shell  structures  is  the  coupling  between  their  bending 
and  stretching  actions.  Through  this  coupling,  a  shell  offers  more  resistance  to  transverse 
loadings  than,  say,  a  flat  plate  of  the  same  dimensions  and  with  the  same  material  properties. 
However,  accrued  to  this  structural  advantage  is  the  price  of  mathematical  complexity.  Even 
in  a  linear  theory,  it  is  often  necessary  to  introduce  additional  assumptions  based  on  plausible 
arguments,  if  the  statics  problem  of  thin  elastic  shells  is  to  lend  itself  to  manageable  solutions. 

The  membrane  theory  of  shells,  for  instance,  assumes  that,  away  from  its  edge(s),  a 

shell  prefers  to  carry  the  applied  loads  by  the  development  of  the  in-plane  stress  resultants 

rather  than  by  the  transverse-shear  and  moment  resultants;  hence,  an  approximate  interior 

solution  can  be  obtained  by  neglecting  terms  associated  with  the  bending  action  in  the  differential 

equations  of  equilibrium.  Such  a  procedure  does,  in  fact,  lead  to  a  good  approximation  of  the 

interior  stresses  for  several  types  of  shells  of  revolution  under  rotationally  symmetric  loads, 

and  loads  varying  sinusoidally  in  the  circumferential  direction  with  period  Ztt  (Refs.  1-3).  On  the 

other  hand,  the  correspondence  between  the  interior  stress  state  and  the  applied  loads  recently 

4  5 

established  by  E.  Reissner  '  for  shallow,  spherical  shells  shows  that  this  is  not  always  true 
for  self-equilibrating  loads.  More  recently,  the  same  correspondence  for  shells  of  revolution 
was  discussed  by  C.  R.  Steele. 

The  present  work  pursues  further  this  subject  of  membrane  vs  inextensional  bending 
interior  stress  state  and  examines  a  larger  class  of  boundary  value  problems  associated  with 
small  deformations  of  shallow,  spherical,  thin,  elastic  shells  (which  is  the  same  as  shallow, 
thin,  elastic,  paraboloidal  shells  of  revolution)  than  those  explicitly  investigated  hitherto.  My 
analysis  differs  from  the  earlier  writers  in  that  the  boundary  value  problenis  will  now  be  solved 
exactly  to  obtain  the  constants  of  integration  associated  with  the  interior  solution  in  terms  of 
the  geometrical  and  material  properties  of  the  shell  and  in  terms  of  the  prescribed  loads  and/or 
constraints.  These  constants  will  then  be  expanded  in  inverse  powers  of  a  large  parameter  to 
unveil  the  nature  of  the  interior  stress  state. 

5  7-9 

Some  known  results  for  small  deformations  of  shallow,  spherical  shells  '  needed  in  the 
subsequent  development  are  recapitulated  in  Sec.  II.  Sections  111,  IV,  V  and  VI  deal  with  shells 
without  surface  loads.  The  shells  are  subjected  to  edge  loads  and  edge  moments  in  such  a  way 
that  the  over-all  static  equilibrium  of  the  shell  is  maintained.  In  Sec.  Ill,  shells  with  prescribed 
edge  loads  and  edge  moments  are  considered.  Once  the  asymptotic  interior  solution  is 
established,  a  comparison  of  the  corresponding  membrane  and  inextensional  bending  stresses 
shows  how  the  nature  of  the  interior  stresses  may  vary  with  the  prescribed  quantities.  The 
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exact  and  asymptotic  influence  coefficients  are  then  obtained.  To  do  this,  we  need  the  explicit 

expressions  for  the  constants  of  integration  associated  with  the  edge  zone  solution.  Since  we 

have  the  explicit  expressions  for  all  the  constants  of  integration,  the  stress  boundary  value 

problem  is  now  completely  (and  exactly)  solved.  As  a  side  result,  it  is  shown  that  the  membrane 

and  inextensional  bending  solutions,  obtained  by  way  of  a  set  of  contracted  stress  boundary 

10 

conditions  established  by  M.  W.  Johnson  and  E.  Reissner  for  shallow,  spherical  shells,  and 
1 1 

later  by  E.  Reissner  for  general  shells,  are  the  leading  term  of  the  asymptotic  interior 
solution.  Shells  with  prescribed  edge  deformations  are  studied  in  Sec.  IV.  Section  V  deals 
with  shells  whose  tangential  edge  displacements  and  transverse  edge  load  and  moment  are 
prescribed.  Shells  with  a  different  type  of  mixed  boundary  conditions  are  treated  in  Sec.  VI. 
There,  the  tangential  edge  loads  and  transverse  edge  deflection  and  moment  of  the  shells  are 
prescribed.  We  then  turn  to  shells  with  surface  loads.  In  Sec.  VII,  shells  subjected  to  axial 
surface  loads  in  the  form  of  polar  harmonics  are  considered.  The  nature  of  the  interior  stress 
state  for  shells  with  various  types  of  edge  support  are  investigated. 

II.  FORMULATION  OF  PROBLEM 

A.  Differential  Equations  in  Polar  Coordinates 


The  system  of  differential  equations  governing  the  small  deformations  of  an  isotropic, 

shallow,  spherical  shell  with  constant  wall  thickness  and  negligible  transverse-shear  deforma- 
7-9 

bility  is 


DV^V^w-lv^F=P^+f 


AV^V^F  +  4  =  -A(l  -  1/) 

K 


(II-l) 


where 


w  =  the  transverse  component  of  the  middle  surface  displacement, 
F  =  a  stress  function  representing  the  direct  stress  resultants, 

R  =  the  radius  of  the  spherical  middle  surface, 

=  the  transverse  component  of  the  surface  load  intensity  vector. 


2,  .  _  3^  .  i  ^  .  J_  9^ 
3r"  >■  '»■  r"  ae" 


■  •  ’,rr  *  I  ",  r  *  T  *  ee 
'  '  r 


i/A  =  Egh,  the  stretching  stiffness  factor, 

D  =  E.  h^/lE(l  —  the  bending  stiffness  factor. 


u  =  Poisson's  ratio, 
h  =  shell  thickness, 

r,  0  =  polar  coordinates  in  the  plane  tangent  to  the  apex  of  the  shell. 


and  where  it  is  assumed  that  the  meridional  and  circumferential  components  P^  and  P^  of  the 
surface  load  intensity  vector  are  derivable  from  a  load  potential  in  the  form 


P 


r 


ar 


i  ^ 

r  30 


(II-2) 


Some  degree  of  nonhomogeneity  is  included  in  the  above  formulation  by  allowing  an  independent 
choice  of  bending  and  stretching  stiffness  factors.  If  the  shell  is  completely  homogeneous,  then 
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=  E,  where  E  is  Young's  modulus.  The  geometrical  properties  of  the  shell  are  shown 

in  Fig.  1. 

The  relevant  stress  resultants  and  couples  are  given  in  terms  of  F  and  w  by  the  following 
relations  (Fig.  2): 


N=-F  +-4F__  +  n 
r  r  ,  r  2  , 00 

r 


=  F  + 
0  ,  rr 


Q  =  —  D(V  w) 


r0  r  ,  0  ,  r 


Qe=-f  (v'w) 


dIw  +t'(—  w  +  w  Q  Q 1 1  , 

[  ,rr  \r  ,r  ^2  ,00/] 

(7  ^  r  +  ^  00)1  > 


M 


=  —  D  I  P' w  + 

■  rr 


M  ^  =  “D{1— i^)  {—  w 

r0  r  ,  0  ,  r 


(11-3) 


The  meridional  and  circumferential  middle  surface  displacement  components  u  and  v  are 
related  to  F  and  w  as  follows: 


u 


r 


w 

R 


r 


,  0 


7u-f  =a(Nq-.n^) 


r 


U  0  +  r(~  v)^  ^  =  2(1  +  p)  AN^q 


(II-4) 


B.  Boundary  Conditions 


To  complete  the  description  of  the  problem,  we  must  supplement  the  differential  equations 
with  an  appropriate  set  of  boundary  conditions.  Throughout  the  present  investigation,  the  shell 
is  to  be  closed  at  the  apex  and  is  to  extend  over  the  region  0^r<:a,  O,$:0^  Zn.  At  the  apex 
r  =  0,  the  shell  is  to  have  finite  stresses  and  displacements.  At  the  edge  r  =  a,  four  consistent 
and  independent  conditions  must  be  prescribed.  Among  the  possible  combinations  of  edge 
conditions  are  the  stress  conditions 


N  =  N  N  ^ 
r  r,  r0 


N  ^  M  =  M  R 
rO,  r  r,  r 


R 


r. 


(11-5) 


where 


.  UlVl  ^ 

R  =  Q  +  -  (ii- 

r  ^r  r  60 

and  N^,  M^,  and  R^  are  the  applied  edge  loads  and  edge  moment;  and  the  displacement 

conditions 


w  =  w  ,  u  =  u  ,  v  =  v  ,  ^  r  ~ 


(11-7) 


where  w,  u,  v,  and  are  the  prescribed  linear  and  angular  displacements  at  the  edge.  In  the 
subsequent  development,  other  possible  combinations  are  also  considered. 
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Fig.  1  .  Geometry  of  middle  surfoce. 
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C.  Solution  to  Homogeneous  Differential  Equations 

5 

The  solution  of  the  homogeneous  system  (II- 1)  is  given  by 

WH=<i'  +  X  ,  Fy=(?;  +  RDV^x  ,  (II-8) 

where  (p  and  ^  are  harmonic  functions,  i.e., 

=  0  ,  (II-9) 

and  X  satisfies  the  equation 

V^V^X  +  - ^  X  =  0  •  (II-IO) 

DAR 

Since  the  shell  is  closed  in  the  circumferential  direction,  the  relevant  displacement  components, 
stress  resultants,  and  moment  resultants  must  be  periodic  functions  of  0  with  period  2n7r, 
n  =  0,l,2,3...  .  Retaining  only  that  f^)ortion  of  the  solution  which  leads  to  the  desired  finiteness 
at  the  apex,  we  have 


oo 

(p  -  Yj 
n=0 


=  y.  r^(B  cos  nG  +  K  sin  n0) 
^  Lj  '  n  n 

n=0 


and 


(II-ll) 


oo 

X  =  ).  [C  J  (kr)  +  H  J  (kr)]  [F  cos  n0  +  G  sinnG] 

'^^‘nn  nn  ^^n  n  ^ 

n=0 

where 

i7r/4  -i7r/4 

,  e  ^  -  e  ^ 

k  -  — ,  k  =  — 

^DAR^  ^DAR^ 


(11-12) 


and  where  A..  B.,  C.,  D.,  E.,  F.,  G.,  H.,  and  K.  are  arbitrary  constants.  J  (x)  is  the  Bessel 
li  i'  i'  i'  1  n 

function  of  the  first  kind  and  of  order  n.  Properties  of  these  functions  can  be  found  in  Refs.  12 
and  13. 

For  the  present  investigation,  it  suffices  to  consider  only  a  special  case  of  (II-ll);  namely. 


(7?  =  A  r^  cos  n0  ,  =  B  r^  cos  n0  , 

^  n  n 

X  =  [H  J  (kr)  +  n  J  (kr)]  cos  n0  ,  (11-13) 

A  I  n  n  n  n  ^ 

where  A  and  B  are  real  constants;  H  and  11  are  complex  constants  with  as  the  complex 
n  n  '  n  n  ^  n 

conjugate  of  to  ensure  real  solutions.  In  the  subsequent  development,  we  shall  replace  r 
by  a  dimensionless  variable  p  defined  by 


r 


a 


(11-14) 
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so  that  0  p  ^  1  and  take  (11-13)  in  the  form  of 

2  n  n 

<P  -  ^  AnP^  cosnG  ,  ip  =  _^i;)  cos  nG 

X  =  [C_ber  (\p)  +  D  bei  (Ap)]  cosnG 

n  n  n  n 


where 


(11-15) 


X  =  -  (11-16) 

\/dar^ 

her  and  bei  are  the  n^^  order  Thomson  functions  and  C  and  D  are  real  constants, 
n  n  n  n 

For  shallow  shells  with  sufficiently  small  bending -to-stretching  stiffness  ratio  DA,  the 
dimensionless  number  X  is  large  compared  with  unity.  For  an  isotropic  homogeneous  shell  of 
constant  wall  thickness. 


A  = 


a 


1  2(  1  —  p'  ) 


(11-17) 


and  A  oo  as  h  0.  In  the  present  work,  we  shall  be  concerned  particularly  with  shells  for 
which  A  »  1.  It  is  known  (Ref.  5)  that  for  this  range  of  A,  the  effect  of  x  is  confined  to  a  narrow 
zone  near  the  edge  of  the  shell;  it  is  therefore  referred  to  as  the  edge-zone  solution.  Away 
from  the  edge,  (p  and  ip  become  dominant;  together,  they  are  referred  to  as  the  interior  solution. 

Corresponding  to  (11-15),  we  have  [cf.  (11-3),  (11-4),  and  (11-6)] 

F  =  I  a^A  p^ - ^  [C  bei  (Ap)  —  D  ber  (Ap)  ]  [  cos  nG 

>  nG 

cos  nG 


a^B  .  1 

brr-^l  [C^ber^(Xp)  +  D^beyxp)]j  cc 


8w 

ar 


naB 


D(1  -  r) 


n  +  ^2^  icherHXp)  +  DheiHXp)] 


a-/b 


n-1 

na(  1  +  v)  AA  p  +  777^7  ;  7m  \ 

n^  RD(n  +  1)  (1  —  i^) 


n+ 1 


+  -  ^  ^  +  1^)  n/~A  ’(Ap)  ] }  cos  nG 

a  ^  n  n  n  n  M 


3^ 

/  .  .  ...  n-1  .  ^  ^n  n+ 1 

na(l  +  ^)AA^p  +  RD(n  +  1)  (1  -V)  P 


n(l  +  u)  nTa 


ap 


C  bei  (Ap)  —  D  ber  (Ap) 
n  n  n  n 


sin  nG 


=  -  [C^bei^(Ap)  -D^ber^(Ap)]  cosnG 

aR  nTA 


= - — -  [C  bei  (Ap)  -  D  ber  (Ap)]  sinnG 

®  R'/A(ap)  "  "  " 
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cos  n0 


where 


N  = 
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n<„  -  1)  Ay-^  t 


"e- 
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re 


M  = 
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M^  = 


n(n-  1)  --  ^ 

\lA  a  1 

n(n-  1)  '  [C^fg^(Xp)  -  D^fg^(Xp)]l  sinnG 


n/A”  (a  p) 


—  n(n  —  1)  B  p  + 

R  nTa 


ICngrc*^^'  “ 


n(n  —  1)  B  p^  ^  +  - - —  [C  g_  (Xp)  —  D  g„  .(Xp)] 

^  n^Gc  ^  n^ed  ^  ^ 


cos  n0 
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re 


n(n-  1)  B  p"'^  +  [C  g  (Xp)  +  D  g  (Xp)] 

R  -v/A  (Xp)  " 


sin  nG 


n  (n  -  1)  B 
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f  (x) 
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n'  X 


bei^(x) 


^0c<^' 


fed<^> 


f  (x) 
sc 


^sd<^> 


f  (x) 
nc 
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n  X  n 
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ber'  (x) 
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2 

^  [berj^(x)  -  i  ber^(x)] 

X 

—  r  “4  beyx)] 
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g  (x) 
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;  bei^(x)  + 


1  - 


2 

ber^(x)  -  ^  ber^(x) 


grd(^> 


gec(x) 


r  2 

=  ber  (x) - -  bei'  (x)  —  —  bei  (x) 

n  X  I  n  X  n 

r  ^ 

=  rbei  (x) - -  ber' (x)  —  —  ber  (x) 

n  X  n  X  n 
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ged<^)  ■  ‘^ber^(x)  + 


i  —  U 


bei’  (x) - bei  (x) 

n  X  n 


g  (x)  =  her’  (x) - ber  (x) 

^sc  n  X  n 


g  ,(x)  =  bei‘(x) - bei  (x) 

^sd  n  X  n 


III.  PROBLEM  OF  PRESCRIBED  EDGE  LOADS  AND  EDGE  MOMENTS 
A.  Prescribed  Conditions  at  Edge 

In  this  section,  we  consider  a  shell  without  surface  loads  acted  upon  by  the  following  system 
of  edge  loads  and  moments  at  p  =  1: 


N  =  N  cos  n0 
r  n 


N  ^  =  S  sinnG 
r0  n 


M  =  M  cos  n0 
r  n 


R  =  R  cos  n0 
r  n 


(III-l) 


where  the  fixed  integer  n  is  greater  than  unity.  A  suitable  solution  to  this  problem  is  that  given 
by  (II~18).  The  boundary  conditions  (II-5)  become 


n(n  —  1)  A  + 

n  2 


[C  f  (\)  -  D  f  ,(\)]  =  -N 
^  n  rc  n  rd  ^  n 


n(n  —  1)  A 


a  n/A 

^ —  [C  f  (X)  -  D  f  ,(X)|  =  S 
2  ^  n  sc  n  sd  '  n 


-n(n-l)B  +  — ^  [C  g  (\)  -  D  g  .(\)]  =  M 
Or  ^  n^rc  n^rd  '  n 


n^(n-  1)  B  +  — - —  [C  f  (X)  -  D  f  ,(X)]  =  aR 
Or  ^  n  nc  '  n  nd  ^  n 


(III-2) 


Solving  (Hi -2)  for  A^  and  we  get 


A  =  — r;  [S  X,  -  N  X,  +  - -  (aR  +  nM  )  X,] 

n  2n(n  —  1)  n  1  n  2  cva  n  n  3 


B 


— ,-^1 -  [(S  +  N  )  X.  +  -^  ^  (aR  X,  -  nM  X,)] 

n  ^  2,  2  .  n  n  4  cva  '  n  1  n  2'  * 

2n  (n  -  1) 


(III-3) 

(III-4) 


where 


X 


R 


1  A, 


2  A, 


1  - 


1  - 


1  —  r 


\ 

n(l  -  L'} 


n(n  +1)  , 

“z-^X—  “3  + 


7^1 


n  +  1  ,  n 

“z-  “3  +-Z  "4 

A 


(III-5) 


X,  = 


3  A, 
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(III-6) 


with 


^4  = 


1 

1  ^ 

r  2n^ 

2 

,  n 

n^(n^  -  1)  (1  -  i^)  1 1 

^  X 

[“2  X 

^  ^  ^  7  ®  4 

^  X 

11 

(n  +  1)  ( 1  —  p') 

Ix 


A^=  1- 
ai  =  -^  (i  =  2,  3,4) 


/  2n  .  \ 

(^2-  —  ^3^-7  ^4) 


(III-7) 


(111-8) 


(3^  =  ber^(X)  bei^(X)  —  bei^(X)  ber^(A) 


^2  =  [ber;^(X)]^  +  [be^  (X) 


=  ber^(X)  ber^(X)  +  bei^(X)  bei|^(X) 


-  ber^(X)  +  bei^(X) 
^4  n  n 


(III-9) 


It  is  clear  from  Fig.  1  that  a  is  the  slope  of  the  middle  surface  at  the  edge  of  the  shell  relative 
to  the  base  plane.  By  the  definition  of  a  shallow  shell  (Ref.  7),  a  «  1  so  that 


1+0? 


(III-IO) 


B.  Some  Relevant  Asymptotic  Relations 

For  large  values  of  X,  we  have  the  following  asymptotic  expansions  for  ber^(X),  bei^(X), 

ber'(X),  andbeiMX)  (Ref.  12): 
n  '  n 


ber  (X) 
n 

|cos 

(— 

TT 

nrr^ 

i  (4n^  -  1)  /X 

+ 

n7r\ 

+ 

0 

n/  ZttX 

8 

00 

f 

c 

0 

X  ^ 

bei^(X) 

1 

sin  1 

X 

TT 

*¥) 

(4n^  -  1)  /  X  _ 

—  + 

n7r\ 

-(?)] 

n/  27rX  1 

■•/2 

8 

8 

2 ) 

ber^(X) 

1 

n/  27rX  1 

cos  1 

(  X 

^-8 

*¥) 

4n^  +  3  /X  V 

nTT^ 

bei’  (X) 
n 

n/  27rX 

|sin 

(— 

I 

(III-ll) 


From  these,  we  get 
a  2 


9 


(III-12) 


a^~  42 


(?)1  ■ 


1  +  o 


We  shall  have  occasions  to  use  these  results  in  the  subsequent  development.  A  new  parameter 
P  is  defined  at  this  point  to  take  the  place  of  A  in  all  subsequent  asymptotic  considerations. 


^  = 


(111-13) 


C.  Asymptotic  Interior  Solution 

For  p  »  n  >  1,  the  following  asymptotic  expressions  for  the  XJs  may  be  obtained  with  the 


help  of  (111-12): 

(n  -  1)  (1  -  i^) 


X,  -  1  + 
1 


X,  ~  1  +  ^  +  o 

3  2p 


(7^ 

(7) 


^2'  ' 


(n  -  1)  (1  -  p>) 

2p 


X .  -  1  +  ^ +  O 
4  Zp 


(7) 

(7) 


(111-14) 


Correspondingly,  we  have  for  p  »  n  >  1 


1 


n  2n(n  ~  1) 
n  —  1 


1  4 


(n  -  1)  (1  -  i^) 


+ 


a  a 


1  + 


2p 

(n  +  1)  ( 1  —  i^) 

2p 


+  O 


(7)1 


1  -  +  o 

2p 


(7)1 


n  2,  2 

2n  (n  —  1) 


1+  (n-l)(l--^)  ^  o 

2p 


(aR  +  nM  ) 
n  n 


(7)]  + 


(111-15) 


X  II.  1)(1  ^Co(^)]  |,_(n  i)(Lji),o^i_jj  I  (III. 


16) 


If  only  the  leading  terms  of  the  expansions  for  the  X^'s  are  retained,  (111-15)  and  (111-16)  are 
further  reduced  to 


\  ~  ^r-r^ — TT  [(S  -  N  )  +  - - -  (aR  +  nM  )] 

n  2n(n  ~  1)  ^  n  n  o;  a  n  n 

R - -  ((S  +  N  )  +  (aR  -nM  ) 

n  n  aa  n  n 


n  2,2  . 

2n  (n  ~  1) 


(111-17) 

(111-18) 


Equations  (111-17)  and  (111-18)  arc  exactly  those  obtained  by  Reissner  for  the  same  problem 
(the  discrepancy  in  the  sign  associated  with  a  is  due  to  choice  of  the  direction  of  the  positive 
middle  surface  normal).  Equations  (111-15)  and  (111-16)  show  that  (111-17)  and  (111-18)  as  well 
as  their  consequences  are  significant  only  if 


and 


(S  +  N  )  +  i^-Li  (aR  +  nM  )  I  «  I  (S  -  N  )  +  - - -  (aR  +  nM  )  | 

n  n  aa.  n  n'  'n  n  a  a  n  n' 


(S  +  N  )  +  IItLJ.  (aR  +  nM  )  I  «  I  (S  +  N  )  +  (aR  -  nM„) 

n  n  aan  n'  'n  n  aa  n  n 


(111-19) 


(III-20) 
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D.  Interior  Membrane  and  Inextensional  Bending  Stresses 


To  examine  the  relative  importance  of  the  membrane  and  inextensional  bending  stresses  in 
the  interior  of  the  shell,  it  suffices  to  consider  two  representative  quantities  and  given 

by 


n<n-l)  .  n-2  _ 


rR 


6M 


h 


6n{n  “  1)  -D 

— ^ — y- — -  B  p  cos  n0 
1  2  n 


From  these,  we  get 
/hA 

,  =o( 

B 


/hA  \ 

—  =0 


(111-21) 


For  p  »  n  >  1,  we  have 
hA 


B 


_n  _  n(n  +  1)  h 


1  + 


—  I’  t  IS„  ^  N„> 


n  -  1 


|1‘  *  «(7)1  ^  °(^)) 


+  - ‘  (1  +  0{-)  I  (aR  +  nM  ) 

aa  '  n  n  n' 


X  (aR^  -  nM  )  +  - iLU - !i]  |i  +  o(-)  ]  (aR  +  nM  ) 

n  n  Z[x  ‘  p  ^  n  n 


{111-22) 


For  a  given  set  of  N  ,  S  ,  R  ,  and  M  ,  the  relative  order  of  magnitude  of  cTt-,  and  cr,.  can  be 
^  n  n  n  n  ^  B  D 

obtained  by  way  of  (111-21)  and  (111-22),  In  what  follows,  we  shall  restrict  our  attention  to  an 
isotropic  and  homogeneous  shell  so  that  [cf.  (111-5)  and  (11-17)] 


a  a  _  a 
h  “  Rh 


n/3(1  - 


(111-23) 


For  the  class  of  problems  for  which  R  =  M  =  O,  (111-22)  becomes 
^  n  n 


hA 


B 


n  _  n( 


n  +  1)  \/3(  1  —  u^) 


— '  [1  +  a-)]  [s  +  N  ] 

*  p  ^  ^  n  n^ 


[1  +  0(-)]  [S  +  N  ] 
‘  p  ^  ^  n  n^ 


(111-24) 


Clearly,  the  nature  of  the  interior  stress  state  depends  on  the  relative  magnitude  of  (S^  +  N^) 
and  (S^  —  N^).  We  now  consider  the  various  relative  magnitudes  and  their  consequences, 

(1)  If  |Sn-  Nnl  «  |Sn  +  Nn|/|x,  we  may  write  (except  for  higher-order 
terms)  (111-24)  as 


hA 

1" 


n  _  n( 


n^  —  1)  (1  —  v)  Jzii  —  u 


(111-25) 


2p 
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Therewith, 


The  membrane  stresses  are  therefore  small  compared  with  the  bending 
stresses  in  the  interior  of  the  shell. 

(2)  If  (Sn  “  Nnl  and  |Sn  +  Nn|/p  are  of  the  same  order  of  magnitude,  both 
terms  in  the  numerator  of  (111-24)  are  equally  important.  Nevertheless, 
we  may  still  conclude  that 


which  is  the  same  as  (1). 

(3)  If  (Sn  —  Nn)  and  (Sn  +  Nn)  are  of  the  same  order  of  magnitude,  we  may 
write  (III-24)  as 


hA 
_ n 

B 


n 


n{n  +  1) 


- 


) 


N 


+  N 


) 


or 


Although  the  bending  stresses  still  dominate  in  the  interior  of  the  shell, 
the  relative  order  of  magnitude  between  membrane  and  bending  stresses 
is  not  the  same  as  that  of  (1)  and  (2). 

(4)  If  +  Nj^)  =  0([Sj^  —  Nj^]/p),  (III-‘28)  continues  to  hold.  But  now  we  have, 
instead  of  (111-29), 


—  =  O(-) 


(5)  If  (S  +  N  )  =  0([S 
n  n  ^  n 


(III-26) 


(III-27) 


(III-28) 


(III-29) 


(III-30) 


(7 

(7 


D 

B 


0(1) 


The  interior  of  the  shell  is  therefore  in  a  mixed  stress  state;  that  is,  both 
the  membrane  and  bending  stresses  are  equally  important. 

(6)  P’inally,  if  (S  +  N  )  «  (S  -  N 
n  n  n  n 


(III-31) 


—  »  1  (III-32) 

The  interior  stresses  are  predominantly  membrane  stresses.  In  particular, 

if  Sn  +  Nn  =  0  (keeping  in  mind  that  Rn  =  Mn  =  0),  there  would  be  no 

inextensional  bending  action  in  the  shell  interior  [see  (III-3)  and  (III-4)]. 

On  the  other  hand,  it  should  be  noted  that  it  is  not  possible  to  have  no 

membrane  action  for  R  =  =  0. 

n  n 

For  the  class  of  problems  for  which  S  =  N  =0,  (III-22)  becomes 
^  n  n 
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(III-33) 


hA 


B 


n  _  n(n  —  1)  "J i{i  ~  v^) 


[1  +0(-)]  [aR^  +nM^] 


— [1  +  0(-)]  (aR  +  nM  ) 
2ui  ^  '  u  '  n  n' 


(1)  If  I  aRj^  +  nMn  I /fi  «  |aRn  — nMn|,  (III-33)  reduces,  except  for  higher- 
order  terms,  to 


hA 


B 


I  Z 

n  _  n(n  —  1)^3(1— r  ) 


aR 
_ n 

aR 

n 


+  nM 

n 

-  nM 


n 


Therewith, 


'B 


In  particular,  if  aR^^  +  ~  ^  (keeping  in  mind  that  =  0), 

there  would  be  no  membrane  action  in  the  interior  of  the  shell  [see 
(III-3)  and  (III-4)  ]. 

(Z)  If  (aRj^  +  nMy^)/p  and  (aRj^  —  nMj^)  are  of  the  same  order  of  magnitude, 
both  terms  in  the  denominator  are  equally  important.  If  their  sum  is 
again  of  the  same  order  of  magnitude  as  each  of  the  terms,  then 


=  0(1) 


Otherwise,  we  must  investigate  the  higher-order  terms  in  the  asymptotic 
expansions  for  Xj  and  to  determine  whether  the  interior  of  the  shell 
is  in  a  mixed  stress  state  or  a  membrane  stress  state. 

(3)  If  |aRj^  — nMnI  «  |  aRn  +  nM^  | /|>jl,  except  for  higher-order  terms, 
we  have 

^^n  _  Zn(n  —  1)  J \  —  r 


Therewith, 


(III-34) 


(III-35) 


(III-36) 


(III-37) 


1 

—  =  O(-)  .  (III-38) 

^B 

Thus,  for  the  class  of  problems  in  which  the  contribution  of  and  to  A^  and  are  negligible, 
the  interior  of  the  shell  is  almost  always  in  a  state  of  inextensional  bending.  The  only  exception 
occurs  when  (aR  +  nM  )/u.  and  (aR  —  nM  )  are  of  the  same  order  of  magnitude  and 

|(aR  -  nM  )  +  (n  -  D  (1  -  v)  (^r  +  nM  )  |  «  I  aR  -  nM  I  .  (III-39) 

'  n  n  Zp  n  n''n  n' 

For  this  exception,  the  relevant  third-order  terms  of  the  asymptotic  expression  for  B^  must 
be  considered  in  order  to  determine  whether  the  interior  is  in  a  mixed  state  or  a  membrane 
state  of  stress.  In  short,  we  may  say  that,  while  self-equilibrating  tangential  edge  loads 
alone  may  lead  to  any  one  of  the  three  possible  interior  stress  states,  it  takes  a  rather  special 
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combination  of  transverse  load  and  moment  (without  tangential  edge  loads)  to  induce  an  interior 
state  other  than  an  inextensional  bending  state. 


E.  Membrane  and  Inextensional  Bending  Solution 

If  only  the  membrane  and  inextensional  bending  components  of  the  solution  are  used  to 

determine  the  interior  state  of  the  shell  without  explicit  reference  to  the  edge-zone  solution, 

the  number  of  constants  of  integration  available  for  the  satisfaction  of  prescribed  conditions  at 

the  edge  is  reduced  to  two:  and  Correspondingly,  a  reduction  of  the  total  number  of 

boundary  conditions  to  be  satisfied  is  necessary.  It  has  been  shown  that  the  appropriately 

10  11 

contracted  stress  boundary  conditions  at  the  edge  p  =  1  are  ' 


(N  - 

-  N  )  +  -  (R 

-  R  )  + 

(M  -  M  ) 

=  0 

r 

r  O' 

r  r 

a  a 

r  r 

,  00 

'^re 

-  N 

rO  a 

(R  -R  ] 
r  r 

© 

+  —  (M  ■ 
a  a  r 

-^r>,0 

(III-40) 


where  the  quantities  N  ,  N  M  ,  and  R  must  now  be  taken  as  the  sum  of  the  solutions 
^  r"  r0  r'  r 

obtained  by  a  membrane  consideration  and  those  by  an  inextensional  bending  consideration. 
For  the  present  problem,  (III-40)  leads  to  the  following  expressions  for  the  constants  and 


1 


B  - 
n 


2n(n  “  1) 
a  a 

2n^(n^  —  1) 


[(S. 


N  )  +  - -  (aR  +  nM  )] 

n  aa  n  n 


[(S  +  N  )  +  (aR  -  nM  )] 

^  n  n  aa  n  n 


(III-41) 

(111-42) 


A  comparison  of  (III-41)  and  (111-42)  with  (111-17)  and  (111-18)  shows  that  the  interior  solution 
obtained  in  this  way  is  just  the  leading  term  of  the  asymptotic  expansion  of  the  exact  interior 
solution.  In  view  of  the  discussion  of  Sec.  III-D,  this  simple  interior  solution  may  not  be  an 
adequate  first  approximation  of  the  exact  interior  solution,  if  (111-19)  and  (111-20)  are  not 
satisfied. 

F.  Influence  Coefficients 

The  deformation  at  the  edge  of  the  shell  will  now  be  expressed  in  terms  of  the  applied  edge 
loads  and  edge  moments.  From  (11-18),  we  have,  for  p  =  1  (or  r  =  a), 

a^B 

w  =  ,  +—  [C  ber  (X)  +  D  bei  (X)] 

n  D(1  -  r)  ^  ^  n  n  n  n  ^ 


= 


naB 
_ n 

D(1  -  1/) 


—  [C  ber'  (X)  +  D  bei'  (X)] 
^  ^  n  n  n  n  J 


3 

^  "“n  ( 1  +  i^)  \  ^/A 

-na(l  +  u)  A  A  +  r. - ^  +  - 

n  RD(n  +  1)  (1  —  i^)  a 


(C  bei’  (\)  -  D  ber'  (X)] 
*  n  n  n  n  ’ 


/ .  .  V  »  »  .  ^  ^n  n(l  +  u)  nTa 

nad  +  r)  AA^  +  RD(n  +  i)  (i  .V) - i - 


(111-43) 
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where  and  are  the  three  edge  displacement  components  and  (S^  is  the  meridional  slope 

change  at  the  shell  edge.  In  order  to  get  the  influence  coefficients,  we  need  and  in  addition 
to  and  To  compute  and  D^,  we  observe  that  and  can  be  eliminated  from  (111-2) 

to  yield 


(C^Ing^^lX)  +  Xf^^(X)]  -  D^(ng^^(X)  +  Xf^^(X)]}  =  (nM^  +  aR^) 


(III-44) 


Solving  (III-44)  for  and  D^,  we  get 


C„  =  -  l(aR  +  nM  )  (f  ,  +  nf  ,)  +  aa{S  +  N  )  (f  ,  +  ^  g  ^)  ] 

rt  A  n  n  rd  sd  rt  m  mH  ^ 


°n  =  -  A 


1 

R  nTa 


n  n  nd  X  ^rd 


((aR  +  nM  )  (f  +  nf  )  -  o-alS  +  N  )  (f  +  ■?  g  >1 
n  n  rc  sc  '  n  n'  '  nc  X  ®rc  ' 


(III-45) 


(III-46) 


where  A  ^  was  defined  by  (111-7).  Substituting  (III-3),  (III-4),  (III-45),  and  (III-46)  into  (III-43), 
we  get  the  following  matrix  equation 

'r 


^wR 

^wM 

KwN 

^wS 

^/3R 

uR 

’^uN 

KuS 

vR 

^vM 

M 


N 


where  the  elements  of  the  influence  coefficient  matrix  are 

"wM  •  =  1^1  =  1^1  "vM 


n  +  1 

-  O'  ,  +  — j 

a  3 


7  “01 


(III-47) 


O'  t  I  .  _  1  —  [  _  n(n  +  1)  _n_ 

\vN  ~  '^wS  ”  *^uR  ~  "\R  n+lA^I  X  [^2  X  ^3  ^2  4 


K ....  =  K  ..e  =  K  =  K  . 


K  —  K  =  I  ^  ^  — — 

uN  ■  vS  'n  +  r  A^ 


[,  ,  2n'^(n  +  1)  (1  -  i^'^)  -  n^(n^  -  1)  (1  - 

r  I 

\ 


1—1^  /  *  1 
(^2  X  ""3 


n^(n  -f  1)^  (1  ~ 
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1  - 


n^(n^  —  1)  (i  —  u)^  +  2n^(n  +  1)  (1  —  u^) 


*  X  '  "2 


,  n^(n  +  1)^  (1  -K^)l  ,  2n^(l  -  u) 

?  I  “^7“  ^ 


1  - 


n(n  +1)  (\  —  v) 


,  n  (1  -  p)  _ 

+  - T -  «3 


1  ~ 


n  (n  -)■  1)  {  \  —  V  ) 


with 


t 


(111-48) 


2n  (n  -  1)  (1  -  i^)  D 

For  |JL  »  n  >  1,  we  have  the  following  asymptotic  expressions  for  the  flexibility  matrix  K; 

n 


K  ~  t 


K 


-  K 
a  m 


K 

n  +  1  p 


K 
a  m 


(-)^  K 
'a'  p 


K 


K 


(n  +  1)  a  m  (n  +  1)  a  m 


n  +  1  ^p 


—  !<■ 
(n  +  1)  a  ^m  'n  +  F  p 


K"  /  ^  \  ^ 

(n  +  1)  a  m  n  +  1  p 


V  ' 


where 


K  =  1  + 
P 


K  =  1 
m 


(n-  1)  (1  -  v) 

2p 

(n  -  1)  (1  - 

Zp 


(111-49) 


It  should  be  pointed  out  that  the  last  two  columns  as  well  as  the  last  two  rows  of  the  matrix  above 
are  the  same.  So,  we  must  consider  either  the  higher-order  terms  in  (111-49)  or  return  to 
(111-47)  if  we  want  to  invert  K  to  get  the  stiffness  matrix. 

IV.  PROBLEM  OF  PRESCRIBED  EDGE  DISPLACEMENTS 
A.  Prescribed  Conditions  at  Edge 

In  this  section  we  shall  again  consider  a  shell  without  surface  loads.  The  shell  is  subjected 
to  a  system  of  edge  loads  and  moments  at  p  -  1  to  produce  the  following  edge  deformations: 


V  =  V  sin  n0 
n 


u  =  u  cos  n9 
n 


w  =  w  cos  n0 
n 


p  =  B  cos  n0 


(lV-1) 


where  n  ^Z.  Again  the  solution  to  this  problem  is  given  by  (11-18).  The  boundary  conditions  (11-7) 
become: 

a^B 


^  [C  ber  (A)  +  D  bei  (A)]  =  w 

D(1  —  i^)  hr  ^  n  n  n  n  ^  r 
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na^B 

5rrr7j  +  —  (C„ber;(X)  +  D^bei;^(X)l  = 


na(l  +  t')  A  A  + 


a^B 


n{l  +  u)  ^fA 


n  RD(n  +  1)  (1  —  i/) 


[C^hei^iX)  -  D^ber^(X)]  = 


—  na(  1  +  t')  A  A  + 


a^B 


n  RD(n  +  1)  (1  —  i^) 


+  <?  +  t'j.XN/T  _  j  ^ 

a  ^  n  n  n  n  ^  n 


Solving  for  A^  and  B^,  we  get 


A  = 


-1 


n  2naA(  1 


n 


^  D(n  +  1)  (1  -  1^) 
Za  a^ 


where 


Y  ,  =  -^  {1  + 
1  ^ 


1_  „  ,  n(nt<yit.) 


‘2  X 


= 


,n  n(n  +  !)(!  +  t') 

O'-,  +  —  a  — 


3  1^2  \  "3 


Y.  = 


4  r3  X  “4 


,  n  n(n  +  1)  (1  +  r) 

O'  ^  +  V  O'  — 


^5  A- 


^6  ZXA^  ^“2  X  “3) 


1  /  n  , 


with 


A2  =  1- 


(n  +!)(!  +  t') 

2X 


/  2n  ,  \ 

("2-  X  “3  +  ^  “4j 


The  quantities  and  X  are  defined  as  in  Secs.  (III-A)  and  (II-B),  respectively. 

B.  Asymptotic  Interior  Solution 

For  p  »  n  >  1,  the  following  asymptotic  expressions  for  the  Y.'s  can  be  obtained 
help  of  (III-12)  and  (III-13): 

V, ■I<<* -I  zo(-L) 


(IV-2) 

(IV-3) 

(IV-4) 


(IV-5) 

(IV-6) 

with  the 
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Y  ~  1  -  . +  o 

2  2fJL 


(j) 


3  2ii 


^  ^  l(n-  1)  +  (n  +  1)  (1  +  P)]  +  O 


(2|Ji) 


Y3~  1. 


a) 

V 

^  1  ^  l(2n-  1)  +  (n  +  1)  (1  +  p)]  + 


±4^  *  o(4) 


(n 


^6-  ^ 


+  +  O^i 


— ^  p  ~  zjl  +  1)  -  (n  +  1)  (1  +  p)]  + 


{IV-7) 


Correspondingly,  we  have, 
,  1 


n  2naA(l  +  r) 

2q'{1  +  r)  w 


^  (n  +  1)  (1  +  t’)  , 
2hl 


°{\)\  *  ‘U‘  * "  ^(rJl  "n 


B 


2p 

2a  (1  +  r)  a/3 

(2tJL)^ 

D(  1  —  r)  (n  +  1) 
2a  a^ 

2a(l  +  u) 

~  2p 

2a  (1  +  r)  a/? 


\a  _l  (n  -  1)  +  (n  +  1)  (1  +  i^) 

1  2h^ 

v)l 

1  + 


1  + 


(2n  —  1)  +  (n  +  1)  (1  +  P’) 


+  O 


m\ 


(IV-8) 


(n 


+  1)  (1  +  i^)  ^  /  1  \1 

1,^)1 

(;7)l 


(U  +  V  ) 

'  n  n' 


1  +  (54_u^  +  o 

2p 


(2p)' 


1  + 


(n  +  1)  ( 1  +  i^)  —  (2n  +  1) 

2p 


+  O 


(IV-9) 


If  the  contribution  from  terms  associated  with  w^  and  can  be  neglected  and  if  terms  of  the 
order  l/p  in  the  remaining  expression  for  and  are  also  discarded,  (IV-8)  and  (IV-9)  are 
further  reduced  to 


(u  -  V  ) 

A  ~  - T - 

n  2na(l  +  i^)  A 


B 


0(1  -  t^)  (n  4-  1) 
Zaa^ 


(u  +  V  ) 
n  n 


(IV-10) 

(IV-11) 


respectively.  Equations  (IV-10)  and  (IV-11)  are  exactly  those  obtained  by  Reissner.  From 
(IV-8)  and  (IV-9),  it  is  clear  that  (IV-10)  and  (IV-11)  are  valid  if 


w  [1  + 
n  ^ 


(n  +  1) 

2p  J 


a/3 

n 

2p 


[1  + 


(n  +  1)  (1  +  i^)  —  (2n  +  1)  1 


«  u  +  V 
'  n  n ' 


(IV-12) 
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and 


«  u  -  V 
'  n  n ' 


(IV-13) 


C.  Interior  Membrane  and  Inextensional  Bending  Stresses 

To  examine  the  relative  importance  of  the  membrane  and  inextensional  bending  stresses  in 
the  interior  of  the  shell,  we  consider  the  ratio 

B 


=°(^) 

'  n' 


(IV-14) 


From  (lV-8)  and  (lV-9),  we  have 
B 


hA 


n(n  +  1)  (1  /  +  y  )  [i  +  0(^)  ]  - 

4/Rh  Inn'  V  ^ 


K-tt) 

X  I,  *  *  o(4)j  -  ^  *  o,i,  -  v„,  |,  X  o(^)] 


2^ 


^  (n  4  IMl  +  P)  +  V  )  [1  +  0(1)  ]  -  £il±Z)  Ifw  -  (i  +  (il 
ZfjL  '  n  n'  ^  V  ^  ^  V  ^  /  [ 


+  1)  i'  +  (n  —  1) 

2]I 


[1  +  o{p 


(IV-IS) 


For  a  given  set  of  u  ,  v  ,  w  ,  and  B  ,  the  relative  order  of  magnitude  of  and  can  be 
^  n"  n  n'  ^  B  D 

determined  by  way  of  (IV-14)  and  (IV-15).  In  what  follows,  we  shall  restrict  ourselves  to 
isotropic  and  homogeneous  shells  so  that 

2 

Hh 


n/3(  1  —  v^) 

Consider  the  class  of  problems  for  which  w  =  B  =  0  so  that  (IV-15)  becomes 
*  n 


B 

_ n 

hA 


=  O 


(u^  +  v^)  [1  +  0(-)] 


(1)  If  |u  +  V  I/p  «  |u  —  V  I,  we  have 
\  f  n'^^  'n  n" 


(IV-16) 


B 

_ n_ 

hA 


Therewith, 


=  O 


/  u  + 

V  \  1 

^  n 

1  2 

U  - 

V  / 

'  n 

n '  J 

^B  1 
—  «  - 


(IV-17) 


(IV-18) 


The  interior  of  the  shell  is  therefore  in  a  membrane  stress  state. 

(2)  If  (uj^  +  Vj^)/p  and  (u^  -  Vn)  are  of  the  same  order  of  magnitude,  unless 
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{IV-19) 


(n  +  1)  (1  +  If) 

2[l 


(u  +  V  )  I  «  I  u  +  V  I 
n  n  '  '  n  n' 


we  have 


<Td 


(IV-20) 


For  the  exception  (if  (IV-l?)  holds),  we  must  consider  the  relevant 
third-order  terms  to  determine  whether  the  interior  stress  state  is  a 
mixed  state  or  an  inextensional  bending  state. 

(3)  If  |u  —V  I  «  |u  +  V  l/yi, 


^  =o<i) 


(IV-21) 


and  again  we  have  a  membrane  interior  stress  state. 

The  above  results  should  be  compared  with  those  given  in  Ref.  5. 

Consider  the  class  of  problems  for  which  u  =  v  =  0.  so  that 

n  n 


B 

_ n_ 

hA 

n 


O 


)] 


(  - 

rn  2|a  ) 


2n  +  (n  +  1)  P’ 


(IV-22) 


It  is  not  difficult  to  see  (with  the  help  of  (lV-22)]  that  the  interior  of  the  shell  for  this  class  of 
problems  is  almost  always  in  a  membrane  stress  state.  The  only  exception  occurs  when 
w^  —  (a^^/2p)  and  a^^/2p  are  of  the  same  order  of  magnitude  while 


(IV-23) 


(2^i) 


For  this  exception,  it  may  be  necessary  to  ascertain  the  relevant  third-order  terms  in  the  various 
asymptotic  expansion  in  order  to  determine  the  interior  stress  state. 

V.  PROBLEM  OF  PRESCRIBED  TANGENTIAL  EDGE  DISPLACEMENTS 
AND  TRANSVERSE  EDGE  LOAD  AND  MOMENT 

A.  Prescribed  Conditions  at  Edge 

Consider  now  a  shell  without  surface  loads  but  with  the  following  constraints  at  the  edge 


V  =  V  sin  nG 
n 


u  =  u  cos  nG 
n 


M  =  M  cos  nG 
r  n 


R  =  R  cos  nG 
r  n 


(V-l) 


where  n  >2.  Substituting  into  (V-l)  the  appropriate  expressions  for  u,  v,  M^,  and  R^  from 
(11-18)  evaluated  at  p  =  1,  we  get 


a^B 


na(l  +  p')  AA  + 


n(l  +  nTa" 


h  RD(n  +  1)  (1  —  p-) 


[C  bei  (X)  -  D  ber  (A))  =  v 
^  n  n  n  n  ^  n 
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(1  +  u)  XsTA 


—  na(  1  +  P')  AA  + 


a^B 


■n  +  ^nVlMfr,)  +  - ^ -  [Cnbei;(X)  -  D^ber;(X)l  = 


-  n(n  -  1)  ^  lCngj.c<^’  “  i  ^ 


"  R  \/A 


n"(n-  1)  [C^f^^(X)  -  D^f^^(X))  =  aR^ 


(V-2) 


Solving  (V“2)  for  A^  and  we  get 


A  =  —  {nM  Z,  +  aR  Z_  +  -rrr - T  fv  Z^  -  u  Z  1 } 

n  2nna  ^  n  1  n  2  A(1  +r)  ^  n  3  n  4^-' 


B 


•  *’  *  >”“„Z5  -  *  MfriT)  l"n  +  "„l  h' 


2\ 


(V-3) 

(V-4) 


where 


with 


7  -  _ L 

1  A. 


1  - 


n^(n^  -  1)  (1  -  1^^)  I  n(l  - 


,  n  —  1  n 

2*^  X  ^3  2^4 

\ 


2  A3 


1  - 


n^(n^  -  1)  (1  -  u^)  I  1  -  u 


■^{h- 


2n^(n^  -  1)  (1  -  1^)  I  i  -  u 


n(n  —  1)  n 

^2  X  ^3  ^4 

A 


n^(n  +  1)  (n^  —  1)  (1  —  v^) 


^2  " 


2n 

X  ^3 


1  - 


n(n  +  1)  (n  -  1)  (1  -  i^)  ]  n 


2X 


+  - o”  O'  - 


•t(i' 


2n^(n^  —  1)  (1  —  i^)  I  1  —  r 


2n  I 

^2“  — ^3  - 


n(n  +  1)  (n^  —  1)  (1  — 

2x'^ 


+  ^0^411- 

X 


n^(n  -i-  1)  (n^  —  1)  (1  — 


X 

^  1  n(l-r)r  (n+1)  ^^n  ,| 

^5  P - ^  ["^2 - ^  «3  ^  X  "4II 


6  A, 


n(n  —  1)  n 

“2 - X —  “s-ir  “4 


^7  =  ^ 


n^(n^  -  1)  (1  -  1^)^ 


1  \  —  V  I  Zr?  \ 

J  X  (“2  X  “3  ■'■  2  “4/ 


A3  =  |1  - 


2n^(n^  —  1)  (1—1^)  I  i  —  V 


- 


n^(n^  —  1)  (n  +  1)  (1  — 

2X^^ 


¥“3  H- 


n^(n^  -  1)  (n  +  1)  (1  -  u^) 

2X^ 


(V-5) 
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B.  Asymptotic  Interior  Solution 


For  fJL  »  n  >  1,  the  following  asymptotic  expressions  for  the  Z.’s  may  be  obtained  with  the 
help  of  (III-12)  and  (III-13): 


Z  ~  1  +  ^  +  o 

1  Z\i 


Z  1  _  ^(n  +  ^  Q 

2,  . 


(-V) 

V  ' 

B) 


1  + 


4n‘^(n  +  1)^  (n-  1)  (1  -  (1  +  u)  ^  qMN 


1  ~ 


1  - 


4n^(n  +  1)^  (n  —  1)  (1  —  (1  +  v) 


(2|jl)- 


Z^~  1 


2(n  —  1)  (1  —  1^) 

(2p) 

2n(n  -H  1)  (1  —  i^) 

(2p)^ 


+  O 


+  O 


1  +  4n^(n^-  1)  (1  +  o 


N 

V 

(^) 
B) 


(V-6) 


At  least  two  terms  have  been  retained  from  the  expansions,  regardless  of  the  order  of  the 
second  term.  It  was  shown  earlier  that  these  second  terms  may  become  significant  under 
certain  circumstances.  Correspondingly  we  have 


A  ~ 


1 


n  ^naa  |  n 


aR  1 


2n(n  +1)  (1  —  i^) 

(2p)^ 


4-  nM  [1  + 
n  ^ 


2(n  -HI)  (1  —  i^) 

2p 


A(1  +  u) 


1  - 


4n^(n  +  1)^  (n  —  1)  (1  —  (1  +  i^) 

(2fi)^ 


A(1  +  u') 


n _  1^  ^  4n^(n  +  1)^  (n  —  1)  (1  —  (1  +  i^) 

(2k)^ 


(V-7) 


B 


-2(n  +  1)  (1  -  ) 

(2fi)^ 


aR  I  1 

n  " 


2n(n  +  1)  ( 1  — 

(2p)^ 


-  nM  [1  - 
n  '■ 


i  A  \  O'  (u  +  V  ) 
2(n  —  1)  ( 1  —  i^) ,  n  n 


2p 


A(  1  +  i^) 


1  + 


A  ^  A  \  I  A  2. 

4n  (n  —  1)  ( 1  “  ) 


.  (V-8) 


We  'shall  also  include  for  later  reference  the  expressions  for  A^  and  retaining  only  the 
leading  term  from  the  expansions  for  the  Z.'s, 


1 


n  2nQ'a 


aiu  —  V  ) 
(aR^.nMn)--^-C^,--;^ 


B 


-2(n  +  1)  (i  -  ) 

(2fi)"^ 


a(u  +  V  ) 
n  n 


(aRn  -  nMj  -  - 


(V-9) 

(V-lO) 


Clearly,  (V-9)  and  (V-lO)  are  valid  only  if 
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n(n  +  1)  ( 1  —  p) 


M  -  ^  f  (u  +  V  ) 

n  \i  a  '  n  n' 


« 


a(u  ~  V  ) 

(all  +  nM  )  -  ■  ,^  ^  P 

n  n  A(  1  +  p-) 


(V-11) 


and 


n(  1  —  u) 

2[i 


(n  -  1)  M  - 
n 


(n  +  1)  aR 


n 


« 


(u  +  V  ) 

(aR  -  nM  )  -  ” 

n  n  A  (1  +  i^) 


Zn(n  -  1)  cv  /  .  ^ 

A  <%  +  '^n> 

(2|jl) 


C.  Direct  Derivation  of  Interior  Solution 


(V-12) 


Although  the  particular  problem  under  consideration  has  not  been  solved  explicitly,  the 
procedure  described  in  Ref.  4  is  clearly  applicable.  In  this  section,  an  approximate  interior 
solution  will  be  obtained  by  the  aforementioned  procedure  to  see  its  relation  to  the  exact  solution 
given  in  Sec.  V-A. 

The  meridional  and  circumferential  displacement  components  u  and  v  can  be  expressed  in 
terms  of  the  stress  function  F  and  the  normal  displacement  w  by  way  of  the  strain-deformation 
relations  (11-4)  as 

u  =  ^  ^  ^dr  -  A(l  +  v)  F  ^  -  (1  +  i-)  RADC  ^ 


(1  +  iv)  RAD?  y 

r  ' 

where 

The  quantities  u^  and  v^,  which  are  the  portion  of  u  and  v  corresponding  to  the  interior  state, 
are  given  by 


H, 


e  2 


A(1  +  i^)  F 


G  i 

—  =  V 


» dr  —  A(  1  +  r)  cp 


Corresponding  expressions  for  M  and  R  in  terms  of  <^,  ip,  x  3^nd  ^  are 


A(1  +  v)  (p 


M  =  ~  D  I  w  +  u 
r  I  .  rr 


M  -  U 
r 


(r  ".r*;?  »,ee)] 
t-(i  -  i-)  (7  »_r  *  ^  x_oe)] 


R 


•D  ((V^w)  +  - - -  (—  w  ] 

^  ,  r  r  r  ,  0  ,  r-^ 


R  1  _  D  (f  +  - - -  (-  X  J  1 

r  ‘  ,  r  r  r  ^ ,  0  ,  r  ^ 


(V-13) 
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where 


M 


i 

r 


ee 


)i 


and 


Rr  =  e’,r*  (V-14) 

are  the  portion  of  and  corresponding  to  the  interior  state.  We  now  approximate  and 
R^  by  discarding  all  terms  involving  x  explicitly  and  we  are  left  with 


.  (V-15) 

(For  a  justification  of  this  approximation,  the  reader  is  referred  to  Refs.  4  and  9.)  At  the  edge 
p  =  1,  we  have  then 

-  ( 1  +  u)  HADt  ^  =  u  cos  n0  ,  v^  -  ^  RAD  ^  ^  ^ 

,  r  n  a  ,  0  n  ' 

M  ^  =  M  cos  n0  ,  R  ^  =  R  cos  n0  .  ( V-1 6) 

r  n  '  r  ,  r  n  ' 

Now  ^  can  be  eliminated  from  these  relations  and  we  are  left  with 

RA(  1  +  P')  R  ^  =  [RA(  1  +  P-)  R  —  u  I  cos  n0 

r  ^  n  n^ 


and 


RA(1  -K  u) 
a 


RA(1  +  i>)  nM 

n 

a 


sin  n0 


(V-17) 


Substituting  into  the  above  two  relations  the  interior  portion  of  the  solution  listed  in  (11-18)  and 
upon  solving  for  A^  and  we  get 


A  =  o 
n  Zna 


—  [(aR 
lo  a  I  n 


aiu  —  V  ) 

+  nM  )  -  ^  ^ 

n  A(1  + 


B  = 
n 


-2(n  +  1)  (1 


V  ) 


(2p) 


1  - 


n^(n^  —  1)  (1  “ 

(2p)^ 


(aR  —  nM  )  — 
n  n 


a(u  +  V  ) 
n  n 

A(1  +  v) 


(V-18) 

(V-19) 


Suppressing  terms  of  the  order  l/p  in  the  presence  of  unity,  (V-18)  and  (V-l9)  are  exactly  the 
same  as  (V-9)  and  (V-10). 


D.  Interior  Membrane  and  Inextensional  Bending  Stresses 


Although  a  detailed  discussion  similar  to  that  given  in  Sec.  IV -C  can  be  carried  out  for 
homogeneous  and  isotropic  shells,  we  shall  limit  ourselves  in  what  follows  to  two  special  cases 
for  the  purpose  of  comparison  with  the  earlier  problems. 

For  a  shell  which  is  free  from  transverse  edge  load  and  moment  so  that  R^  =  M^  =  0,  we 
have,  from  (V-7)  and  (V-8),  for  (p  »  n  >  1), 


(V-20) 
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where 


K  =  n^(n  +  1)^  (n  -  1)  (1  -  u)^  (1  +  v) 


Thus,  similar  to  shells  with  only  prescribed  tangential  edge  loads  {Sec.  III-D),  the  interior  of 

the  shell  may  be  in  a  membrane,  a  mixed  or  an  inextensional  bending  stress  state  depending  on 

the  relative  magnitude  of  the  quantities  (u^  +  v^)  and  (u^  —  v^). 

On  the  other  hand,  if  the  shell  is  constrained  in  the  tangential  directions  so  that  u  =  v  =0 

n  n  " 

we  have  (for  p  »  n  >  1) 


(aR^  -  nMJ  +  [(n  -  1)  - 

(aR  +  nM  )  +  t.  ~  '''>  (m 
n  n  ^  \ 


n  +  1 


aR 


aR  V 

_ n\ 

n  2[i  ) 


(V-21) 


Thus,  similar  to  shells  with  prescribed  transverse  displacement  and  meridional  change  in  slope 
at  the  edge  (along  with  vanishing  tangential  edge  displacements),  the  interior  of  the  shell  is 
almost  always  in  a  membrane  state  (Sec.  IV-C).  The  only  exception  occurs  when 


aR  +  nM 
n  n 


(V-22) 


and 


|(aR  +  nM  )  +  +  1)  <1 - !l]  /jyi  _  |  «  |  aR  -  nM  |  ,  (V-23) 

'  n  n  |JL  \n2}j./''n  n' 

in  which  case  the  relevant  third-order  terms  in  may  have  to  be  considered  in  order  to 
determine  whether  the  interior  is  in  a  mixed  or  an  inextensional  bending  stress  state. 

Note  that  these  results  are  in  agreement  with  our  earlier  observation  (Sec.  III-D).  They 
demonstrate  once  again  that  while  self-equilibrating  tangential  edge  loads  alone  may  lead  to 
any  one  of  the  three  possible  interior  stress  states,  it  takes  a  rather  special  combination  of 
self-equilibrating  transverse  edge  load  and  edge  moment  to  induce  an  interior  stress  state 
other  than  an  inextensional  bending  state. 

VI.  PROBLEM  OF  PRESCRIBED  TANGENTIAL  EDGE  LOADS 
AND  TRANSVERSE  EDGE  MOMENT  AND  DEFLECTION 

A.  Prescribed  Conditions  at  Edge 

In  this  section,  we  consider  a  shell  with  the  following  edge  conditions  at  p  =  1: 


N 

=  N 

cos  n0 

N  ^  =  S  sin  n©  , 

r 

n 

r©  n 

M 

=  M 

cos  n©  , 

w  =  w  cos  n©  , 

(VI-l) 

r  n  n 


where  n  >2.  Note  that  (VI-1)  differs  from  (III-l)  only  in  the  last  condition;  the  prescription  of 

the  transverse  force  R^  in  (III-l)  is  now  replaced  by  that  of  the  transverse  deflections  w.  For 

fixed  values  of  N  .  M  ,  and  S  ,  the  results  of  this  section  and  those  of  Sec.  Ill  are  identical,  if 
n  n  n 

the  R^  is  taken  to  be  exactly  that  amount  of  transverse  force  needed  to  produce  the  transverse 
displacement  w^.  On  the  other  hand,  the  explicit  asymptotic  interior  solution  for  (VI-1)  will 
make  some  rather  interesting  conclusions  easily  accessible. 
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Observing  (11-18),  (VI-1)  becomes 


- r  B  +  — ^  [C  ber  (\)  +  D  bei  (X)]  =  w 

D(1  —  r)  n  ^  n  n  n  n  '  n 


-n(n  -  1)  B  +  - —  [C  g  (A)  -  D  g  A\)]  =  M 

n  n^rc'  n^rd'  n 


n(n  -  1)  A  +  —  [C  f  (X)  -  D  f  ,(X)]  =  -N 

n  'v/A  ^  n  rd  ^  n 


n(n  -  1)  A - [C  f  (X)  -  D  f  .(X)]  =  S 

n  n/A  ^  sc  n  sd  ^  n 


Solving  (VI-2)  for  A^  and  B^,  we  get 


A  =  — 
n  aa 


[m^  +  n<n-l)^n  -»-)  D 


B  = 


1  -  u 


2  )  T  -S  T  ] 

n(n‘^  -  1)  " 

(n  +  1)  Dw 


n  n  +  1 


T .  +  M  -  aa(S  +  N  )  T, 
4  n  5  '  n  n'  6 


where 


with 


B.  Asymptotic  Interior  Solution 

For  |ji  »  n  >  1,  the  following  asymptotic  expressions  for  the  T^'s  may  be  obtained 
help  of  (III-12)  and  (III-13): 


(VI-2) 


(VI-3) 

(VI-4) 


(VI-5) 


with  the 


26 


-p  1  +  - !i) 

M  2^ 


T  ~  1  -  ^  .  -  (  ^  +  o 

iZ^i) 


T,  ~  1  +  2n(n^l)  (l^y)  ^  Q 

U\^) 

T,  ~  1  +  ^<n-  1)  (1  -£)  +  o 
T.  ~  ^  [l-^  [(2n_l),2,,  +  o(-^)] 


(7) 

(7) 

b) 


^  (2m)^ 


‘  (2„)  = 


1  —  V 

1  +  +  O 

Z\1 


fe)l 


Correspondingly,  we  have 


A - ^ 

n  aa 


“n  *  1-  *  0(;,)|  * 


X  ,„N„  |.  -  .  o(A,)j  _  s„  (.  2  *  °(> 


and 


B  -  (1  -  r) 
n 


2n(n  —  1 )  (1  —  u) 


ZU 


1  +  +0 

1  -  .  ^/ 1 


(2p) 

4aa 


(7)] 


2n  -  (1  -  Zu) 
Z\x 


+  O 


5 


(n  +  1)  (Zfi)- 


1  ^  ^°[-^)\  (Sn  +  N,)j 


If  only  the  leading  term  of  the  expansions  is  retained,  (VI-7)  and  (VI-8)  become 


A  - - — 

n  cv  a 


n(n-  1)  (1  -  .)  D  „  +  _g 

?  nj  n  n'j 


and 


B  ~  (1  -  12)  —  w  +  ^  y'*  M - ^^-5 - ^  (s  +  N  ) 

"  "  (2^)2  "  (n+l)(2^)^ 

respectively.  Clearly  (VI-9)  and  (VI-10)  are  valid  as  long  as 


/  A  \  /  A  \r\  (S+N) 

^  n(n-  1)  (1  -  .)  D  ^  n  - n 

n  2  n  n  +  1  2u. 

a 


« 


M  +  n(n-  1)  (1  ->-)  D  w 
n  2  n 


(nN  -  S  ) 
n  n 


n(n  -  1) 

n(n  —  l)(l-'r)D  2n— (1  —  21^) 


w  — 
n 


2p 


M  -  - —  (S  +  N  ) 

"  2^^(n+l)  "  " 


« 


D  ,,  ,  iin±^  M  - 

an  /o  \2  n 


4cva 


(2p) 


(n  +  1)  (2p) 


-5  (S  +  N  ) 
i  n  n 


(VI-6) 


(VI-7) 


)1 

(VI-8) 


(VI-9) 


(VI-10) 


(VI-11) 


(VI-12) 
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C.  Interior  Membrane  and  Inextensional  Bending  Stresses 

From  (VI-7)  and  (VI-8),  we  have 


_  Qfa(l  —  r) 
hA  ""  h 


X  IVI  - 
n 


/Dw  M  V 

r  /  A  \ 

1 

2"+  2 

1  +  o(-^) 

'  a^  2|ji^' 

1  V  ' 

J 

n(n  —  1)  (1  —  r)  2n  —  (1  —  2r) 

2  2  .3 

a  4|jl 


^  (S  +  N  )  [1  +  O(-) 
3  n  n  ‘  Li 


(n  +  1)  2|jl 


2m. 

n(n  -  1)  (1  -  r)  D 


w  +  M 
n  n 


X  [1  +  O(^)  ]  + 


-  1) 
1 ,  ,1 


InN  -  S  1 
^  n 


1  +  O 


(A)i  -  ff- 


-  i') 

2|jL"(n  +  1) 


X  (N  +  S  )  [1  +  0( 
n  n  ^  |JL 


(VI-13) 


A  rather  unique  feature  of  the  present  problem  now  reveals  itself.  In  all  the  problems  considered 
earlier,  if  we  set  all  but  one  of  the  prescribed  quantities  equal  to  zero,  the  interior  stress  state 
remains  the  same  for  each  problem  regardless  of  which  quantity  is  nonzero.  This  is  not  so 
for  the  present  problem. 


If  M  =  N  =  S 
n  n  n 


0. 


B 

n 

hA 

n 


n( 


n  —  1)  J  —  V  ^ 


(VI-14) 


whence  in  the  interior  of  the  shell. 

If  w  ^  S  =  N  =0, 
n  n  n 


B 

_ n_ 

hA 

n 


1  -  V 


n/3(1  -v" 


[1  +  o(:7)i 


(VI-15) 


ag  and  are  of  the  same  order  of  magnitude  in  this  case. 


D 

If  w  =  IVI  =  S  =0, 
n  n  n  ' 


^  ,  (n-  1)  (1  ^  0(1, 

hA  I - -=-  '  V 

"  «i3(l  - 


{VI-16) 


Therewith 

a 


=o(i) 


(VI-17) 


and  we  have  a  membrane  interior  stress  state. 

In  arriving  at  (VI-14),  (VI-15),  and  (VI-16),  we  have  assumed  that  the  shell  is  homogeneous 
and  isotropic  so  that 

2 


aa. 

h 


n/3(1  -u‘ 


(VI-18) 
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Under  the  same  assumption,  we  now  consider  the  class  of  problems  for  which  w  =  M  =0 

n  n 

For  these  problems,  (VI-13)  reduces  to 


B 

_ n_ 

hA 

n 


n(n  —  1)  (1  —  v) 


X 


(nN  -  S  ) 
'  n  n 


(S  +  N  )  [1  +  0(-)l 

_  n  n  ^ ^  p  ^ _ 

'  +  o(^)]  -  (S^  +  NJ  [1  +  0(1)] 


(VI-19) 


All  three  stress  states  are  possible  in  the  interior  of  the  shell,  depending  on  the  magnitude  of 
(Sn  +  ^  ^n^‘  Equation  (VI-19)  should  be  compared  with  the  corresponding  result  given 

in  Ref.  4. 


VII.  SHELLS  WITH  POLAR  HARMONIC  AXIAL  SURFACE  LOADS 
A.  Solution  to  Differential  Equations 

Let  the  surface  load  intensity  vector  be  axial  with  magnitude 
P^(r,  e)  =PQp'^cosne  (n^l) 


Within  the  context  of  shallow  shell  theory,  an  axial  and  a  normal  surface  load  are  completely 

5 

equivalent.  Thus  we  may  write 


=  P^  =  ^0^^  cos  n0  (n  >  1) 

P  s  =  0  .  (Vll-1) 

r  0 

The  particular  solutions  to  (Il-l)  with  P^  as  given  in  (VIl-1)  are 


F 

P 


w 

P 


^0^  n+2 

“  T7 — FTT  P  cos  n0 
4(n  +  1)  ^ 

PqAR^P^  cos  n0 


(Vn-2) 


Although  the  particular  solution  for  w  has  the  same  form  as  the  inextensional  bending  component 
of  the  homogeneous  solution,  and  therefore  may  be  omitted  from  further  consideration,  it  will 
be  kept  throughout  the  subsequent  development  to  facilitate  later  comparisons  with  existing 
results. 

The  corresponding  particular  solutions  for  the  resultants,  couples,  tangential  displacement 
components,  and  meridional  slope  change  are; 

Npp  =  cosne  ,  Nqp  =  PqRp"  cos  nG 

N„  =  -  X  PnRp”  sinnG  ,  Q  =Q  =0 
r0p  4  0^  *  ^rp  0p 

2 

M  =-M„  = -n(n  -  1)  (1  -  i/)  P^p^"^cosn0 

rp  0p  ^  0 
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=  n(n  -  1)  ( 1  -  1^)  — 


R^p  =  n^(n  —  1)  ( 1  —  i')  ^ 


n(  1  +  v) 


P^aARp^^^  sinnG  , 


p  4(n  +  1)  0 

(n+2)(l  +  i/)^  ^^n+1  ^ 

""p  =  ”“4(ri  +  1) -  ' 


^  PqAR^p'^'^  cosne  (VII-3) 

The  complete  solutions  to  the  problem  can  now  be  obtained  by  summing  these  particular  solutions 
and  the  corresponding  homogeneous  solutions  given  in  (11-18).  The  four  constants  of  integration 
which  appear  in  the  resulting  expressions  are  to  be  determined  by  the  constraints  at  the  edge 
p  =  1.  For  n  =  1,  our  results  are  just  those  obtained  by  Reissner.^'^  Observe  that  in  this  case, 
the  particular  solutions  for  the  transverse  shears  and  the  bending  moments  vanish  identically. 
Moreover,  the  interior  solutions  for  these  quantities  and  for  the  in-plane  stress  resultants  also 
vanish  identically. 

In  the  next  few  sections,  the  nature  of  the  interior  solution  for  shells  with  various  types  of 
support  will  be  considered.  In  particular,  we  shall  show  that,  to  the  first  approximation,  the 
interior  stresses  of  a  shell  with  a  clamped  edge  and  one  with  a  simply  supported  edge  are  the 
same,  and  that  the  results  from  a  membrane  and  inextensional  bending  analysis  are  equivalent 
to  this  first  approximation  of  the  interior  solution.  For  n  =  1,  the  membrane  solution  turns  out 
to  be  the  exact  solution  to  (II-l)  for  a  shell  fixed  tangentially  but  free  otherwise.  The  interior 
of  the  shell  in  all  cases  is  dominated  by  the  membrane  stresses.  On  the  other  hand,  for  n  >  2, 
the  interior  of  a  shell  with  free  edge  is  primarily  in  a  state  of  inextensional  bending. 


9w 

P  _ 

~w  ■ 


B.  Clamped  Edge 

Consider  first  a  shell  with  a  clamped  edge  so  that 

u  =  V  =  w  =  =  0  (at  p  =  1)  .  (VII-4) 


Substituting  into  (VII-4)  the  appropriate  expressions  for  the  four  displacement  quantities,  we 
get 


a^B 


D(1  -  1^) 


|C  ber  (X)  +  D  bei  (X)]  =  —  P„AR 
'  n  n  n  n  0 


na^B 
_ n 

U{1  -  i') 


+  (C  ber'  (X)  +  D  bei'  (X)]  =  -nP.AR^ 
^  '  n  n  n  n  '  0 


a^Bn 

na(l  +  r)  AA^^  +  +  i)  (i  _ 


n(  1  +  i^) 


[C  bei  (X)  —  D  ber  (X)] 
^  n  n  n  n  ^ 


n(l  v) 
4(n  +  l) 


P^ARa 


*  The  particular  salution  given  in  Ref.  5  far  the  radial  displacement  component  is  actually  the  expression  for  the 
meridional  displacement  component. 
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X(1  +  t')  n/X 


a^B 


—  na(l  +  v)  AA  + 


n  RD(n  +  1)  (1  —  v) 


iC^beV^iX)  -  D^ber|^(X)] 


(n  +2)  (1  +  u) 
4(n  +  1) 


Solving  (VII-5)  for  A  and  B  ,  we  have 
n  n 


Y^  +  nY3-8(n+l)(Y3-nY^)l 

P„a^(l  -  u^) 

- TTi—  |Y4-4(n+l)(Y^-nY^)]  , 

(2fi) 

where  the  Y.'s  are  just  those  given  by  (IV-5).  If  |jl  »  n  >  1,  we  have 


A 


n  4n 


1  -  +  o 

Zil 


b)l 


B 


P^a^d  - 


1  -  <”  +  --1}  +  O 

d\l 


(^)l 


If  terms  of  the  order  i/\i  are  also  suppressed,  (VII-8)  and  (VII-9)  become 


^0^  1 
A  ~  (1  +  a-)] 

n  4n  ^  V 


B..  ~  - 5 -  (1  +  0(^)1 

(2kl)^ 


n 


For  p  »  n  >  1,  we  have 
2 


ap,  ^-’VRh  Aj 


2  2 

If  the  shell  is  isotropic  and  homogeneous,  a  /Rh  =  0(p  ),  and  therefore 


—  ■  odl 


Thus  the  interior  of  the  shell  is  in  a  membrane  stress  state. 


C.  Simply  Supported  Edge 

If  the  shell  is  simply  supported,  the  conditions  (VII-4)  are  replaced  by 
u  =  v  =  w  =  =  0  (at  p  =  1) 


Substituting  into  (VII-14)  the  appropriate  expressions  for  the  displacement  components 
bending  moment  for  the  present  problem,  we  have 


D(1  -1.) 


+  — —  (  C  ber  (X)  +  D  bei  (X)]  =  -P„AR 
^  >  n  n'  n  n  '  0 


2 


(VII-5) 

{VII-6) 

(VII-7) 

(VII-8) 

(VII-9) 

(VII-10) 

(VII-11) 

(VII-12) 

(VII-13) 

(VII-14) 
and  the 
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1 


,  n(n  -  1)  (1  -  I')  P  a 

n(n  -  1)  - —  (C  g  (X)  -  D  g  ,(X)]  - - - - 


2 

a  a  B 


na(l  +  i^)  AA  + 


n  D(n  +  1)  { 1  — 


n(l  +  i^)  4A  ^  ^ 

a  ”  4(n  +  1) 


2 

cva  B 


—  na(  1  +  1^)  AA  + 


n  D{n+l){l— i^) 


-Q -  +  VA  [c  bei'  (X)  -  D  bei'(X)] 

a  ^  n  n  n  ^ 


(n  +  2)  ( 1  +  i^)  P^ARa 
4{n  +  1) 


{VII-15) 


Solving  {VII-15)  for  A^  and  B^,  we  get 


A  = 


PqR 


n  8  n(  n  +  1 )  A  , 


nK^  +  (n  +  2)  K2 


4(n  +  1)  ^  4n(n  —  1)  (1  —  i') 


(VII-16) 


P  a^  (1  -v^) 

Bn  =  -  ^TTT- -  (K3) 

(2(jl)  a. 


(VII-17) 


where 


K  =a  -  +  (n  +  1)  (1  +  1^),  (n  +  1)  (1  )  (  _  H  „  ) 


_  1  -  n(n  +  l)(l-i^)/  n^  ^ 

*^2  -  "4  -  - n -  <“3  -  X  °'4> 


K, 


1  _  n^(n^  —  1 )  ( 1  —  I  _  \  —  u 


n(n  —  1)  n 

^3  .2^4 

A 


K  =  Q,  +  R  a  -  BIP  +  1)  (1  +  »^) 

4  “3  X  4  ,2 


Kc  = 

5  4 


1  + 


2n  (n  -t  1)  {  \  —  u) 


-  1  [d  -  i;)  +  2(n  +  D]  +  ^ 


I  1 

la-2  -  ^  0-3] 


^6  =  «4  + 


n^(n  +  1)  (1  -  1-^)1  1 


zx 


^  [(1  - 1-')  + 


(n  +  1)  (1  +  i^) , 


The  cv^’s  are  defined  by  (III-7). 

For  \i  »  n  ^1,  we  have  the  following  asymptotic  expressions  for  A^  and  B^ 


A 


n  4n 


1  - 


(3  *-1^) 
4p 


+  O 


(d)l  ■ 


(VII-18) 
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(VIl-19) 


^  *  °(^)1 

If  terms  of  the  order  l/[i  are  also  neglected  in  comparison  with  unity,  (VII-18)  and  (VlI-19) 
become 

1 

An  ~  (1  +  0<-)]  (VII-20) 

and 


B 


1  - 


(n  +  1 


B 


n 


P^a^d  - 


ii  +  a^)] 


(Vll-21) 


A  comparison  of  (Vll-18)  and  (VII-19)  with  (Vll-8)  and  (VlI-9)  suggests  that,  to  the  first 
approximation,  the  interior  state  of  the  shell  is  the  same  whether  the  edge  of  the  shell  is 
clamped  or  simply  supported.  To  say  it  in  another  way:  for  shells  with  a  small  bending-to- 
stretching  stiffness  ratio  so  that  p  »  n  ^  1,  the  interior  of  the  shell  is  rather  insensitive  to 
the  difference  between  the  two  types  of  edge  supports.  It  is  interesting  that  the  second-order 
correction  terms  in  (VlI-18)  and  (Vll-19)  are  exactly  half  the  corresponding  terms  in  (VlI-8) 
and  (VII -9). 


D.  Membrane  Analysis 

Since  both  types  of  edge  support  considered  in  the  last  two  sections  lead  to  a  membrane 
interior  state,  we  might  suspect  that  the  corresponding  solution  obtained  by  way  of  a  membrane 
analysis  would  provide  a  good  approximation  to  the  exact  solution  for  shells  with  a  clamped  or 
simply  supported  edge.  We  shall  presently  show  that  this  is  so. 

Let  us  consider  now  a  shell  with  vanishing  bending  stiffness.  The  governing  equations  for 
this  membrane  shell  can  be  obtained  from  (ll-l)  by  setting  D  =  0.  Thus  we  have 

V^F  =  -RP  =  -  RP^p"  cos  ne 
m  n  0'^ 


V^w  =-RAV^V^F 
m  m 


(VIl-22) 


The  solutions  to  (VlI-22)  for  and  w^  are  (keeping  in  mind  that  the  shell  is  closed  at  the  apex 
and  that  the  stress  resultants  and  displacements  are  finite  there) 


=  la^A 


P.a^R 
0  n+2 

4(n  +  1)  ^ 


cos  n0 


a  -y  I 

mn  n  ,  ^  n\  ^ 

5(T-7)  P  ^  p  I 


(VII-23) 


Correspondingly,  the  tangential  displacement  components  are 

n+ 1 


/II  ^  A  A  n-1  ^  ^mn^ 

nad  +  d  AA^^p  +  RD(n  +  1)  T~7) 


n(l  +  P')  aP^ARp 


n+1 


4(n  +  1) 


sin  n0 


33 


m 


3t^  n+1 

.  a  B  p 
/ .  .  »  A  A  n-1  . 

—nail  +  u)  AA  p  +  rrFTT — .  . v  / ^ - : 

'  '  mn'^  RD(n  +  1)  (1  —  i^) 


(n  +  Z)  (1  +  v)  aP^ARp 


n+li 


4(n  +  1) 


cos  n0 


(VII-24) 


The  constants  and  are  determined  by  the  usual  membrane  (tangentially)  fixed  edge 

conditions  (Fig.  3). 


(at  p  =  1) 


(VII-25) 


The  result  is 


PqR 

4n 


B 


Pga^d  -vh 


I-7I-I722I 


(VII-26) 


Thus  we  see  that,  except  for  the  forces  and  moments  associated  with  the  inextensional  bending 
of  the  shell  (which  are  identically  zero  since  the  shell  offers  no  resistance  to  bending),  the 
membrane  solution  provides  a  good  first  approximation  to  the  interior  solution.  If  the  portion 
of  the  above  results  associated  with  is  regarded  as  the  inextensional  bending  component 

of  the  complete  solution  so  that  the  transverse  shears  and  bending  moments  can  be  derived  from 
the  expression  for  w  (D  0  idv  this  purpose),  a  comparison  of  (VII-26)  with  (VII-10)  and 
(VII-11)  or  with  (VII-20)  and  (VII-21)  shows  that  the  result  from  such  an  analysis  is,  in  fact, 
the  leading  term  of  the  asymptotic  expansion  of  the  interior  solution  for  shells  with  either  a 
clamped  edge  or  a  simply  supported  edge. 

E.  Shells  Fixed  Tangentially  But  Free  Otherwise 

In  this  section,  we  show  that  a  membrane  interior  state  is  still  ensured  if  the  last  two 
conditions  in  (VII-4)  are  replaced  by  R^  =  =  0.  Moreover,  if  n  =  1,  we  have  the  interesting 

result  that  the  membrane  solution  turns  out  to  be  the  exact  solution  of  (II-l).  Inconclusive 
though  they  may  be,  these  results  provide  at  least  a  partial  and  quantitative  substantiation  of 
the  appropriateness  of  the  edge  support  used  in  the  membrane  theory  of  shells. 

Consider  then  a  shell  with  vanishing  tangential  displacements  and  transverse-shear 
resultant  and  moment  (Fig.  3). 

u  =  V  =  =  R^  =  0  (at  p  =  1)  .  (VII-27) 

For  the  present  problem,  (VII-27)  takes  the  form 
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na(l  +  u)  A  A  + 


a^B 


"n  ■  RD(n  +  1)  (1  -^.)  i  ~  "  4(n  +  l)  ^0 

-na(l  +  .)  AA^  +  RD(n!lMl-.)  +  H ' 


ARa 


(IL_+  ^)  (1  -^)  p  ARa 
4(n  +  1)  ^0  ^ 


P  a 

-n{n  -  1)  B„  +  — !—  (C  g^^{\)  -  D  e  .(X)  ]  =  n{n  -  1)  { 1  -  .;) 


n  ^  n®rc  n®rd' 


P  a^ 

n^{n  -  1)  B„  +  (C^f^^{X)  -  D^f^^{x)  )  =  n^(n  -  1)  (1  -  ,;) 
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n 


R  ^/A 


Solving  (VII-28)  for  and  we  get 


(VII-28) 


A  = 


^0 ^  [ 4n^(n^  —  1 )  (1  —  u) 


n  8n(n  +  1) 


(Z^  +  Z^)  +  nZ^  +  (n  +  2) 


(VII-29) 


P„a^  {I  -i^^)  f  ,  2,  2  . 

B  =-— -  _  gn_(n  -  1)  (1  -u)  .  _z)  +  Z 

4^4  1  ^4  (^3  Z,)4Z, 

where  the  Z^'s  were  given  by  (V-5). 


If  H-  »  n  ^  1,  we  have 


PR 

A - ^  |1  +  O 

n  4n 


i^)\ 


{VII-30) 


{VII-31) 


B 


P^a^  (1  -i/^) 


1  +  O 


(;?)| 


{VII-32) 


{VII-33) 


"  {2m.)  ^  'M 

For  a  homogeneous  and  isotropic  shell  (again  \x  »  n  ^  1) 

u  '  ^j.  '  V 

which  is  the  first  part  of  our  contention.  Note  that  correction  terms  in  (VII-31)  and  (VII-32)  are 
0(1/^"^)  while  these  same  correction  terms  are  0(l/p)  in  Secs.  VII-B  and  VII-C. 

For  n  -  1,  (VII-29)  and  (VII-30)  become  [cf.  (V-5)) 


A,  = 


PqR 


(VII-34) 


and 


B,  = 


-P^a^  (1  -r^) 
(2m)‘‘ 


(VII-35) 
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respectively.  Moreover,  the  last  two  conditions  of  (VII-28)  become 


C,g^^(X)-D^g^^(X)  =  0 


or 


Cl  =  Di  -  0 


(VII-37) 


Equations  (VII-34),  (VII-35),  and  (VII-37)  constitute  the  second  part  of  our  contention.  It  is  not 
difficult  to  verify  that  the  shell  is  still  in  a  state  of  over-ail  static  equilibrium. 

F.  Shells  Fixed  Transversely  But  Free  Otherwise 


As  a  contrast  to  (VII-27),  we  consider  next  a  shell  supported  in  such  a  way  that  (Fig.  4) 


N=No  =  w=  |^=0  (atp  =  l) 
r  r0  9r  \  / 


(VII-38) 


Observing  (11-18),  (VII-38)  becomes 
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n  _  2 


n(n  —  1)  A - y 

n  2 


a  nTA 
nA 


\fA 


(C  f  (X)  -  D  f  .(X)]  =  P„R 
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na  B  .  2 

^  +  —  [C  ber'  (A)  +  D  bei'  (A)]  =  -nP.AR 
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(VII-39) 


If  n  >  2,  we  may  solve  (VII-39)  for  A^  and  B^  to  get 


_  PqR 

n  “  4(n  +  1)  A, 
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(VII-40) 
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B  =  — ^ - 
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(VII-41) 


with 


A  ^n  ,  n 

■^8  ““2  X 


I -71-1733  I 


(VII-42) 
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The  a^'s  are  defined  by  (III-7).  For  |jl  »  n  >  1,  we  have 
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n 


4(n  +  1) 


B 


n 


P^a^d  -  v) 
(n  +  1)  (2|ji)^ 


For  n  =  1,  {VII~39)  becomes 

X  1  1 
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{Vll-45) 


The  first  two  conditions  of  {VIl-45)  are  the  same;  we  are  then  left  with  three  relations  for  the 
three  unknown  B^,  and  appearing  in  these  relations.  The  solution  for  B^  is  again  given 
by  (VII-41)  with  n  =  1. 

With  all  the  prescribed  edge  constraints  satisfied,  we  are  still  left  with  an  arbitrary 
constant  which  appears  only  in  the  expressions  for  the  tangential  displacement  components 
u  and  V.  A  closer  examination,  however,  reveals  that  the  contribution  of  the  terms  associated 
with  A^  to  u  and  v  is  exclusively  in  the  nature  of  rigid  body  motion. 

In  all  cases,  as  long  as  p  »  n  ^  1,  we  have,  for  a  homogeneous  and  isotropic  shell, 


Although  the  membrane  stresses  again  dominate  in  the  interior,  the  relative  importance  of  the 
two  stress  states  as  given  by  {VII-46)  is  more  marginal  than  that  exhibited  by  the  shells 
previously  considered  in  this  section. 


G.  Free  Edge 

Since  several  types  of  edge  constraints  have  all  led  to  a  membrane  interior  state  for  the 
axial  surface  load  under  consideration,  we  would  naturally  ask  whether  there  is  any  set  of  edge 
constraints  which  would  lead  to  an  inextensional  bending  interior  state.  To  answer  this,  let  us 
consider  a  shell  with  a  free  edge  so  that 

N=N,^  =  M=  R=0  (atp  =  l)  .  {VII-47) 

r  rO  r  r  ' 

Clearly,  this  is  a  permissible  set  of  edge  conditions  only  if  the  applied  load  is  self-equilibrating. 
Therefore,  we  must  restrict  the  integer  n  to  greater  than  unity.  Observing  (11-18)  and  (VIl-3), 
(VII-47)  becomes 

n(n  -  1)  A  +  -  (C  f  (X)  -  D  f  ,(X)]  =  P.R 

n  a^  n/A  ^  n  rd'  J  4  0 
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"  R^^A: 
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Solving  (VII-48)  for  and  we  get 
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where  the  X.’s  are  just  those  given  by  {III-5).  For  |ji  »  n  >  1,  we  have 

I*  t  ^  *  o(  _ 
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n  ^  2,  2  ,, 

4n  (n  —  1) 


1  -  +  O 

2n 


(^)l 


If,  in  addition,  the  shell  is  homogeneous  and  isotropic,  we  have 


t  -  °m  -  °(5) 


Thus  the  interior  of  the  shell  is  primarily  in  a  state  of  inextensional  bending. 

Vm.  SUMMARY  AND  REMARKS 


(VII-48) 


{VII-49) 


{VII-50) 


(VII-51) 


{VII-52) 


{VII-53) 


In  the  foregoing,  we  have  solved  a  series  of  boundary  value  problems  exactly  and  investigated 

the  corresponding  asymptotic  behavior  of  the  so-called  interior  solutions  subject  to  the  assump- 

4  5 

tion  \x  »  n.  For  shells  without  surface  loads,  the  present  work  supplemented  known  results,  ' 
with  the  explicit  determination  of  the  second-order  corrections  to  the  leading  term  of  the 
asymptotic  interior  solutions.  These  correction  terms  enabled  us  to  establish  the  conditions 
under  which  the  leading  terms  provide  an  adequate  first  approximation  to  the  exact  interior 
solutions.  They  also  led  us  to  some  refinements  of  the  correspondence  between  the  interior 
stress  state  and  the  boundary  conditions.  Exact  and  asymptotic  influence  coefficients  were 
obtained.  For  shells  with  polar  harmonic  axial  surface  loads,  our  results  showed  that  the  shell 
interior  is  in  a  membrane  state  for  the  various  types  of  edge  support  considered.  On  the 
other  hand,  if  the  edge  is  free,  the  interior  of  the  shell  is  in  a  state  of  inextensional  bending. 

It  is  hoped  that  these  results  may  contribute  to  a  better  understanding  of  the  interplay  between 
the  interior  membrane  and  inextensional  bending  stresses. 
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